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Abstract 

An analytical theory is developed for parametric interactions in metamaterial multilayer struc- 
tures with simultaneous nonlinear electronic and magnetic responses and with near-zero refractive- 
index. We demonstrate theoretically that electromagnetic fields of certain frequencies can be 
parametrically shielded by a nonlinear left-handed material slab, where the permittivity and per- 
meability are both negative. The skin depth is tunable, and even in the absence of material 
absorption, can be much less than the wavelength of the electromagnetic field being shielded. This 
exotic behavior is a consequence of the intricate nonlinear response in the left-handed materials 
and vanishing optical refractive-index at the pump frequency. 

PACS numbers: 78.66.Sq, 42.65.Ky, 68.65.Ac, 42.25.Bs 



The recent advancement of fabrication technologies has opened up numerous avenues 
in designing metamaterials, which are composite structures with modified electromagnetic 
properties. The main thrust of metamaterial design is to push the limits of ever increas- 
ing control over the effective material parameters, permittivity, e and permeability, /i. By 
accessing previously forbidden regions of material parameter-space, many highly unusual 
electromagnetic properties and fascinating phenomena have been predicted to occur, not to 
mention the possibility of novel device applications that may not be realized with conven- 
tional materials 

It is well known that over certain frequencies, typical metals can reflect electromagnetic 
(EM) fields and can thus be used as electromagnetic shielding materials. In response to 
the incident EM field, free electrons in the metal collectively oscillate and generate fields of 
opposite phase relative to the incident field, leading to an exponential decay of the total fields 
inside the metal. These evanescent electric (E) and magnetic {H) fields cannot propagate 
due to the negative sign of e and positive /i, intrinsic to metals at the appropriate bandwidth. 
The skin depth, i.e., the distance EM fields drop to 1/e of their surface value, is determined 
primarily from the frequency of the source, and the corresponding material parameters. If 
the situation arises whereby a structure possess simultaneously negative effective e and /i, 
a type of metamaterial referred to as left-handed material (LHM), the EM fields can then 
propagate inside the material with a negative refraction at the interface, and ideally the 
system exhibits transparency. Thus, conventional linear LHMs cannot be used to shield 
electromagnetic fields. 

In this letter, we demonstrate that this picture is drastically modified when nonlinearity 
of the magnetic response is taken into account, creating a controllable shielding effect in 
LHMs, accompanied by a parametric reflection. Unlike metals, that behave as passive 
reflectors or mirrors, a nonlinear metamaterial slab can function as an active mirror that 
changes the incident field frequency upon reflection. We show that the nonlinear-induced 
skin depth, 6, which is a measure of the shielding strength in LHMs, is tunable and much 
less than the wavelength. A, of the EM fields being shielded, even in the absence of material 
absorption. It is vastly different in conventional metals, where the skin depth is typically 
fixed and on the order of A. This shielding effect is the consequence of backward parametric 
interaction due to the nonlinear magnetic response of the LHM and near- zero refractive index 
at the pump frequency. Recently, many interesting phenomena associated with near-zero 
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refractive index metamaterials have been reported, such as directive emission 

and stephke transmission!^. In the more conventional materials (positive e and yu), also 
called right-handed materials (RHMs), most nonlinear effects originate from the electronic 
response. With left-handed materials however, the magnetic permeability is on equaHooting 
and plays a crucial role. A nonlinear magnetic response can indeed arise in LHMs 
and exact phase matching can be achieved between forward and backward waves 
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12|. In the majority of past works, the 



leading to enhanced second-harmonic generation 
nonlinear response was restricted to be either electric or magnetic. Moreover, to simplify 
the problem, the slowly varying envelope approximation (SVEA) was often used. In this 
paper, without recourse to SVEA, a general theory is developed for parametric interactions 
in metamaterial multilayer structures simultaneously possessing nonlinear electronic and 
magnetic responses. By solving Maxwell's equations exactly, multiple reflections at the 
RHM/LHM interfaces are accounted for, thus properly describing the nonlinear parametric 
interaction between forward and backward waves. 

Consider a ID A^-layer periodic structure with period d = dr + di, where dr and di 
are, respectively, the layer thickness of the RHM and the LHM. The total length of the 
structure is L. The 2;-axis coincides with the direction normal to the plane of the layers. 
The structure is configured such that the right-handed medium has a nonlinear electronic 
response, whereas the left-handed medium has a nonlinear magnetic response. The linear 
permittivity in the RHM, e^-, is given by the standard Drude-Lorentz form, 

e.H = i+ , ""r ■ - (1) 

where ujpr is the effective plasma frequency, Ur is a resonant frequency, and 7 is the damping 
factor of the medium. The linear permeability in the RHM, /i^, is constant. In the LHM, 
the linear permittivity, e/, and permeability, /i/, are similarly described by, 

ei{u) = 1 + — |L^,^^(^) = 1 + (2) 

where Upi is the effective plasma frequency, ui and ujq are the resonant frequencies of the 
medium, and F is a geometrical filling factor. The parametric interaction is a three-wave 
mixing process in a nonlinear medium that involves a strong pump beam and two relatively 
weak beams, the signal and idler. The distinction between signal and idler is arbitrary, just 
as long as energy and momentum are conserved. In general, due to its dispersive nature, 



LHM possess left-handed properties only in certain frequency regimes. We assume the signal 
(ui) is a negative- index wave, while the idler {lj2 < t^i) and the pump {u^ = ui + 002) are 
both positive-index waves. We consider TEM modes propagating along the 2;-direction as 
an x-polarized E field and y-polarized H field. In the RHM, the electric field of the signal, 
El, and idler, E2, are thus described by 



<PE2 ^ul 



(3) 



Similarly, the corresponding magnetic fields in the LHM are given by, 

, 9 H — ^eiz/^ii-ni = —1^11^2 



^ H T^2lfJ'2lH2 = —K21HI , 



(4) 



where Kir = {'^i/c)'^{87^f^irX^e^E3) and nu = {uJi/c)'^{87ceiiXTn H3), for i = 1,2, are nonlinear 

(2) (2) 

coupling parameters that couple the signal and idler waves. The coefficients, Xe and Xm 
are the second-order nonlinear electric and magnetic susceptibilities, respectively. Here E^ 
and H3 represent the pump wave. With the assumption of plane waves and negligible pump 
depletion, Eqs. ([3]) and (jlj) can be solved analytically. In particular, we consider metama- 
terials with near-zero refractive-index at the pump frequency. This approach stems from 
the fact that the pump photons carry negligible momentum in three-wave mixing processes, 
and thus forward and backward coupling between signal and idler become dominated in the 
LHM, leading to enhanced parametric reflection. Although the refractive index is close to 
zero, the pump can still be coupled into the structure through matched impedance methods 

nn 

or tailoring the waveguide appropriately |5|]. 

The fields inside each layer are linear superpositions of forward and backward waves. 
Each forward and backward wave contains two modes due to the parametric interaction 
between the signal and idler. The general EM solutions in the nth layer can be succinctly 
written in terms of the modal fields, ^^^J (z) (for q = r,l): 

Ei''\z) = <!>^;J{z)+C,rKi'\^), 

(5) 



e;'^-\z) = c;M:^{z) + <i><;\z 
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in the RHM. For the LHM, we have 



(6) 



The couphng coefficients Ciq and C2q {q 



r,l) are given by Cig = KiJ{r]g 
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), and 



C2q = K2q/{f3f - ^ig), where = {ui/ cfeir^iir and = {uji/cyeajju, for i = 1,2. The 
modal fields are written in terms of plane waves, 



Inserting Eqs. ([5]) and into the nonlinear wave equation yields the propagation constants 



When the nonlinear parameters Kiq and K2q vanish, the propagation constants (3q and r]q 
correctly reduce to their linear counterparts for the signal and idler, i.e. Pq = kig (sig- 
nal) and rjq = k2q (idler). Furthermore, the solutions in Eqs. ([5]) and (El) reduce to linear 
signal and idler waves. When the nonlinear parameters are non-zero, Eqs. and 
contain all possible parametric interactions among forward and backward waves between 
signal and idler. The modal coefficients a^jjj, a^^^, ^tq/!,! and a^^^ in Eq. ([7]) can be found 
by matching boundary conditions that require the continuity of tangential E and H at 
each interface. The inherent nonlinear boundary conditions can be avoided using Maxwell's 
equations: H = c/{iujfir)dzE in the RHM, and E = c/ {iujei)dzH in the LHM. To pro- 
ceed, we define state vectors in the RHM and LHM: la^''"'*) = ('^r//3' '^Ife^' '^ri!^'*' '^r/'^^)"^' 
and = (f^;/^, O/t)], ft/fe^"*, ct/*)^^)^- The transformation between state vectors occurs via: 

|a,(")) = Tr\ar^), and |a^("+i)) = Tj|a/")), yielding |a^("+i)) = TzT^|a^(")). The transfer 
matrices and Ti are calculated using the Maxwell relations above. After some tedious 




=Vi exp(i/5g2) + aj^^ exp(-z/5gz), 
=a-^^^e^^{iriqz) + aJ^J exp(-ir7gz) . 
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f3q and r]q, 
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where, (pg = exp{ii3qdg), and 6q = exp{ir]*dq) . The generahzed impedances in the right- 
handed and left-handed media at the signal {n = 1) and idler {n = 2) frequencies are 

written as Znr^ = (u;n/inr)/(c/?„r), ^nrr, = (uJnfinr) / {cT]nr) , in RHM, Znip = {c(3nl) / {i^nenl) , 



(g = r, /). Similarly, 



in LHM. Moreover, (5\q 
I 



l3q, I32q = [3*, r]iq = 7]*, and r]2g = r]q, 
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and, 
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The input and output layers are arranged to be RHM, and the signal and 
idler (co'2) waves can be incident from both the right and left sides of the structure, 
as shown in Fig. [H By matching tangential fields at the input and output interfaces, 
the transmission and reflection can be calculated via, \U) = A~^{\Sfi} — M\Si,)), 
where \U) = {Eif^, E2fji, Eibi, £2^1)'^ . The left source, Sl, is written as \Sl) = 



6 



E , E 

l/I 2/1 
— ► 



E , E 



E ,E 

m 2(R 



E ,E 

IbR IhR 



FIG. 1: (Color online) Metamaterial multilayer considered in this paper. The amplitudes of the 
left input fields are represented by -Ei/l and E2fL, and the right ones by Euji and E2bR- The 
amplitudes of the left output fields are represented by -Ei^l and -E26L i and by Ei fji and E2 /r for 
the right side of the stack. The subscripts 1 and 2 refer to the signal and idler waves, respectively. 

{EifL, EifiZ^^ , E2ji^Z2^^)^ , and the right source, Sr, compactly as ISr) = 
{EibR, —Eib^Z^j^, E2i,ji, —E2f,^Z2j^^)^ , where the ambient impedances at the frequencies 
ui and UJ2 for the right and left sides of the structure are, Zi^ = {oJiHiRj / {ckiz) and 
ZiL = (uJifiiL) / (ckiz) . The 4x4 matrix A is given by 

^ -R.M,L,^M.L. (13) 



R2 M3L1 + M4X2 



where, 
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-2iL 0/ \o -z;,' 



(14) 



The 2x2 submatrices Mi {i = 1, 2, 3, 4) arise from the 4x4 matrix M, 

M = P^TI;' = M , (15) 

where Pr and Ir are given by Eqs. ([9]) and (fT2|) . The transfer matrix T expresses the modal 
coefficients in the A^th layer in terms of the first one: la^^-* 

In the following, we consider the pump wavelength to be operating at A3 = 1.2 /xm 
(corresponding to ^3 0), and take F = 0.77, uq = 120 THz, Upi = 300 THz, Upr = 150 THz, 
uJr = 200 THz, uJi = 120 THz, and 7 = 0. We begin by investigating parametric reflection 
and transmission through a single slab of nonlinear LHM. Figure [2] shows the transmission 
and reflection normalized to the incident electric fleld at the idler wavelength, A2 = 3.6 /xm. 
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FIG. 2: (Color online) Transmission and reflection of signal (Ei) and idler (£'2) for a single nonlinear 
LHM slab as a function of slab thickness (left panel) and Xm H^/ fj,2i (right panel). The incident 
field is at the idler wavelength (A2 = 3.6 /um) and the pump is set to 1.2 /im, thus the signal is 
Ai = 1.8^m. At the idler wavelength, ij.21 = 1.72. In the left panel, Xm is fixed at 0.06 

while in the right panel, the thickness of the slab remains at 1 ^m. 

There is no signal frequency at the input; the signal is generated inside the stack due 
to the nonlinear effect of the materials. Thus, the signal reflection [Ei) exhibited in the 
plots implies backward transmission. The frequency of the signal, ui, is determined by 
the frequencies of the idler and pump via energy conservation, uJi = 00^ — 002, yielding 
Ai = 1.8 fim. The left panel shows that the signal and idler reflections increase with slab 
thickness, while transmission decreases. The right panel illustrates the variation of the 
reflection and transmission with the pump field when the thickness of the slab is fixed at 
1 yUm. 

To demonstrate electromagnetic shielding, we now consider a three-layer structure com- 
posed of RHM sandwiched between two layers of LHM. An electromagnetic field at the idler 
wavelength A2 = 3.6 fim is incident from both sides of the stack. Shown in Fig. [3] are the 
intensities of the electric (top panel) and magnetic (middle panel) fields normalized to the 
input electric field. The z component of the Foynting vector inside the stack for the signal 
and idler is also shown in the bottom panel. As before, there is no signal frequency at the 
input; the signal is generated inside the stack where the backward parametric interaction 
between the signal and idler in the negative-index layer leads to an exponential decay from 
the input interfaces. After a slight penetration, the electromagnetic fields are virtually zero. 
Thus, a nonlinear negative-index layer can be used as a parametric mirror that changes fre- 
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FIG. 3: (Color online) Parametric shielding by a trilayer metamaterial stack (see text). The 
thickness of the LHM and RHM layers are 0.5 /im and 1.5 ^m, respectively. The incident field 
at the idler wavelength (3.6 ^m) impinges both sides. Both signal (solid blue) and idler (dashed 
green) decay exponentially inside the stack. The skin depth is ~ 0.1 ;um, with xe = 0.01e2r, 

(2) 

Xm H3 = 0.06/X2Z, where €2r = 1-68 and fi2i = 1.72 at the idler wavelength. 
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FIG. 4: (Color online) The results when Xm = 0. Other parameters and inputs are the same as 
those in Fig. [3j The fields are significant inside the stack. The absence of EM shielding is clearly 
evident. 



quency upon reflection via the nonlinear interaction in the penetration layer. Remarkably, 
the skin depth can be 30 times smaller than the wavelength of the incident field, as seen by 
the approximate expression. 



S 



V2 \/\fIllfJ^2l\ / Af 



|ei//^i/| 



+ 



(16) 



showing explicitly that S is inversely proportional to the pump field strength. The dramatic 
effects of the nonlinear shielding vanish if the nonlinear magnetic response of the LHM is zero, 
as exhibited in Fig. HJ where the same quantities as those in Fig. [3] are plotted except with 
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Xm = 0. The top two panels of Fig. H] show that the electromagnetic fields are significant 
inside the stack. The bottom panel shows the Poynting vector, where the signal generated 
in the nonlinear RHM (the middle section from 0.5 fim to 2 fim) propagates straight through 
the linear LHM (the two end sections) without idler interaction. The small but finite energy 
fiow arises from the nonlinear effect in the RHM. For a purely linear structure, a signal wave 
cannot be generated, and the idler becomes a standing wave with net zero power flow inside 
the stack due to oppositely directed input flelds with identical amplitudes. 
This work is supported by NAVAIR's ILIR program sponsored by ONR. 
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